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We study several issues related to the different choices of time available for the classical and 
quantum treatment of linearly polarized cylindrical gravitational waves. We pay especial attention 
to the time evolution of creation and annihilation operators and the definition of Fock spaces for the 
different choices of time involved. We discuss also the issue of microcausality and the use of field 
commutators to extract information about the causal properties of quantum spacetime. 
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I. INTRODUCTION 

The quantization of polarized gravitational cylindri- 
cal waves has received a lot of attention in recent years 
Blllllllllllilili- This is partly due to the fact that 
this system provides a tractable, yet non-trivial, reduc- 
tion of full general relativity and hence is an ideal frame- 
work to explore several issues involved in the quantiza- 
tion of gravity. Some intriguing phenomena, related to 
the existence of large quantum gravity effects, have been 
discussed by studying precisely this model 0, S IE S 13 ■ 
It has also been argued that some manifestations of quan- 
tum gravity, such as the smearing of light cones are, in- 
deed, present and can be understood in this simplified 
setting 0. 

One of the interesting points behind the obtained re- 
sults is the realization of the fact that the physical Hamil- 
tonian is a function of the free field Hamiltonian for 
a 2-1-1 dimensional, axially symmetric, massless scalar 
field evolving in an auxiliary Minkowski background 
P, 0, 01 . As we show in the first section of the pa- 
per, this free Hamiltonian naturally appears when one 
linearizes the system, thus suggesting that, in a precise 
sense, it can be considered as the free part of an inter- 
acting model. However the full interacting Hamiltonian 
is obtained by adding a very specific type of terms to the 
free part, namely, just functions of it. Here we plan to 
explore the consequences of this functional relation be- 
tween the two physically relevant Hamiltonians for the 
system and explore how this affects the causal structure 
of quantum spacetime. To this end we will pay attention 
to the smearing of the light cones due to quantum gravity 
effects within the framework of linearly polarized cylin- 
drical waves. Instead of considering the full information 
encoded in the metric tensor we will concentrate on the 
causal structure provided by light cones. An interesting, 
albeit somewhat indirect way, to look at this structure is 
to study the commutators of field operators at different 



spacetime points. These are the basic objects to discuss 
the commutativity of observables and the microcausal- 
ity of the model; conventionally, a physical model should 
be such that observables commute for space-like inter- 
vals. This has been discussed in the standard perturba- 
tive quantum field theory framework for simple examples 
such as scalar or fermion fields (see e.g. In fact, 

microcausality is one of the key conditions to prove such 
important results as the spin-statistics theorem 0, . 



Here we will use the commutator of the scalar field that 
describes linearly polarized cylindrical waves as a way to 
get information about the causal structure of quantum 
spacetime. As we show later it is possible to give exact 
expressions for this commutator both for the evolution 
provided by the free and the full physical Hamiltonians. 
We will use these expressions to study in a quantita- 
tive way the smearing of the light cones as a function 
of the three-dimensional gravitational constant and ex- 
plore some physical issues, in particular the appearance 
of singularities as a consequence of having a Hamiltonian 
bounded from above [l4.ll5|. 
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The rest of the paper is structured as follows. In Sec. 
nil we discuss how the free Hamiltonian is derived from 
the linearized cylindrical wave model. Sec. IIIII deals with 
the classical and quantum dynamics of cylindrical gravi- 
tational waves under the evolution provided both by the 
free Hamiltonian and the physical Hamiltonian. We will 
pause here to discuss and compare on a familiar example 
(the harmonic oscillator) the main features of the time 
evolution defined by functionally related Hamiltonians, 
both from the classical and the quantum points of view. 
This will provide valuable insights for the problem con- 
sidered in this work. Sec. IIVI is devoted to the study of 
microcausality. We will look at the main features of the 
field commutators and study the smearing of the light 
cones due to quantum gravitational effects. We end the 
paper with a discussion of the main results and perspec- 
tives for future work. 
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II. CYLINDRICAL WAVES IN LINEARIZED 
GRAVITY 

Linearly polarized cylindrical waves in general relativ- 
ity can be described by the spacetime metric 0, |^ : 



(1) 



where Z e M is the coordinate of the symmetry axis and 



ds| is the three-metric 



-N^dt^ + e<{dR + N^dtf + r^de^ 



(2) 



From this three-dimensional point of view, R € M"*" and 
6 ^ correspond to polar coordinates, is the radial 
component of the shift vector and N is the lapse function. 
All metric functions (7, r, iV, and N-^) depend only on 
the time and radial coordinates, t and R. 

Unless otherwise stated, we adopt a system of units 
such that c = H = 8G3 = 1, where c is the speed of 
light, H is Planck constant, and G3 is the effective Newton 
constant per unit length in the direction of the symmetry 
axis Ql- In these units , th e gravitational action of the 
system has the form |^ 0| 



5 = / dt 
Jti 



dR {p^j + PrT -\- p^ip) — H 



(3) 



where the dot denotes the derivative with respect to t, the 
p's are the momenta canonically conjugate to the metric 
variables and Ti. is the total Hamiltonian 



n = 2 1 



-7oo/2 



poo 

/ dR [NC + N"'C"] . (4) 
Jo 



The first term is a boundary contribution at infinity 
[700 ■= l{R 00)] and the second term is a linear com- 
bination of the Hamiltonian constraint C and the (radial) 
diffeomorphisms constraint C^: 



-7/2 



2r' 



ir' 



PyPr + ^ 



c = 



The prime denotes the derivative with respect to R. The 
gauge freedom associated with these constraints can be 
removedb^ imposing, respectively, the gauge fixing con- 
ditions 



a Dy : 



X 



■.= r-R = 0, 



X :=P7 = 0. 



On the other hand, the Lagrangian form of the action 
can be obtained from the relations between momenta and 
time derivatives of the metric provided by the Hamilton 
equations, 

p^N = -e'''/^r + e'^/^N^r', 
p^N = e''^^ri:~e-''^^N^r^\ 

p^N = + 2 (iV«e-^/2y . 



From a three-dimensional perspective, the system de- 
scribes an axially symmetric model consisting of a scalar 
field V' coupled to gravity 0, the line element being ||2J). 
A particular classical solution is a vanishing scalar field in 
three-dimensional Minkowski spacetime or, equivalently, 
Minkowski spacetime in four dimensions. In this solution, 
N = 1 and r = R, whereas the rest of metric fields and 
momenta vanish (i.e. ip = j = = p^ = p^ = pr = 0). 

In this section, we will consider this solution as a back- 
ground and discuss first-order perturbations around it. 
In other words, we will analyze the linearized theory of 
gravity around this Minkowski spacetime, as it is usually 
done in the perturbative, quantum field theory approach 
to gravity. In order to expand the metric fields around 
the classical solution, let us call 

r = R + f, N = l + N. 

Up to first-order terms in the fields, the expression of the 
three-metric becomes 

dsj = -{l + 2N)dt^ + 2N'^dtdR+{l+-/)dR'^ 
+ {R^ + 2Rf)de'^, 

while the four-dimensional metric is given by 

ds^ (1 - ^j)dsl + (1 + ij)dZ^. (5) 

Here, it is understood that the product of ^ with any 
other metric field vanishes in the perturbative order con- 
sidered. On the other hand, regularity on the axis of 
symmetry imposes the following conditions: 

7(i? = 0) =0, N^{R = 0) = 0, 

f{R = 0) = 0, f'{R = 0) = 0. (6) 

In order to discuss the linearized gravitational system, 
we must keep up to quadratic terms in the fields in the 
action A straightforward calculation leads to the 

result 



S 

n 



dt 



dR 



dR (p^7 + Pf-f + p^^) — Ti. 



P4, , i?(^') 



l\2 



2R 



P^Pf + NC + N"-C 



+ (2-7oo)r'(i?-^(X3). 



Here, Pf '■— Pr is the momentum canonically conjugate 
to f, and the linearized constraints are: 



7 , 



= Pr - 2p' 



The diffeomorphisms gauge freedom can be fixed just 
like in the cylindrical reduction of general relativity, 
namely, by demanding that r = R + f coincides with 
the radial coordinate Q- We thus impose the gauge fix- 
ing condition := f = 0. It is easily checked that the 
Poisson brackets {x'^, C^} of this condition with the lin- 
earized constraint do not vanish, so that the gauge fixing 
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is admissible. Dynamical consistency of the gauge fixing 
procedure requires, in addition, 



R — 



Hence, the shift must be chosen as = pj. Finally, 
the momentum conjugate to r is fixed by solving the 
diffeomorphisms constraint: pf — 2p^y. In this way, the 
canonical pair (f, Pf) is removed from the set of degrees 
of freedom. The action of the resulting reduced model is 



^1 



dt 



dR 



(^Pyi +p^'4' - NCl 



dt [Ho-p^^iR ^ oo) = 0)] 



Here, Ci = —7' is the Hamiltonian constraint of the re- 
duced linearized system, and 



dR 



4 
2R 



(7) 



Remarkably, the condition employed to eliminate the 
Hamiltonian gauge freedom in full cylindrical gravity 
can be used as well to fix the corresponding gauge in 
the linearized theory. The gauge fixing x '■— P-^ — 
is acceptable, because the Poisson brackets of x and Ci 
differ from zero. In addition, consistency of the chosen 
gauge demands the vanishing of 



X 



dRNCi } = -N' 



Therefore, N has to be independent of the radial coordi- 
nate. Actually, we can set TV = by demanding that the 
total lapse equals the unity at spatial infinity. On the 
other hand, taking into account the regularity condition 
©, the solution to the Hamiltonian constraint is simply 
7 = 0. This allows us to remove the canonical pair (7,^7) 
from the system and arrive at a constraint-free model in 
linearized gravity. 

The degrees of freedom of this system are the field ip 
and its momentum. The reduced action is 



^2 



dt 



-Ho 



dR p^ip 



Note that Hq, given in Eq. {T)), is the Hamiltonian of 
a massless scalar field with axial symmetry in three- 
dimensions. Furthermore, in the gauge that we have 
selected, the three-dimensional metric of the linearized 
gravitational theory is exactly that of Minkowski space- 
time and contains no physical degrees of freedom. The 
scalar field ip determines the norm of the Killing vector 
dz, and appears in the four-dimensional metric of the 
gauge-fixed, linearized model in the form |(SJl, but with 
(is| substituted by the flat metric 



-dT^ + dR^ + R^dO^ 



(8) 



where we have renamed T the time coordinate of the 
reduced system. 

Summarizing, the perturbative description provided in 
linearized gravity for cylindrical waves with linear polar- 
ization around four-dimensional Minkowski spacetime is 
equivalent to a massless scalar field with axial symmetry 
in a three-dimensional flat background. The dynamics of 
this field is dictated by the free Hamiltonian i7o, which 
generates the evolution in the Minkowskian time T. 

It is worth noticing that the action and metric of the 
gauge-fixed model in linearized gravity reproduce in fact 
the results that one would obtain from the gauge-fixed 
model in full cylindrical gravity by working just in the 
first perturbative order, i.e., by keeping in the action and 
metric, respectively, up to quadratic and linear terms in 
the field and its momentum. In this sense, the gauge 
fixing and linearization procedures commute. 



III. TIME COORDINATES AND EVOLUTION 
FOR CYLINDRICAL WAVES 



Systems with Functionally Related 
Hamiltonians 



One of the significant features of gravitational cylin- 
drical waves is the existence of two distinct, physically 
relevant Hamiltonians (or equivalently two distinct time 
coordinates) to define both the classical and quantum 
evolution. One is the Hamiltonian Hq [given in Eq. iQ] 
that generates the dynamics in the linearized gravita- 
tional theory; the other is the Hamiltonian H that pro- 
vides the energy per unit length along the symmetry axis 
in general relativity llj,|l5|. In fact, they are func- 
tionally dependent, since H = 2(1 — e~^°^^). In order 
to gain insight into the relation that can be established 
between the dynamics associated with these two differ- 
ent Hamiltonians, we open this section by discussing a 
similar situation in a very simple example provided by 
the harmonic oscillator. 

The usual description of the harmonic oscillator in a 
phase space coordinatized by {xq,Po) comes from its stan- 
dard Hamiltonian hQ{xQ,po) — {pq + uj'^Xq)/2. The dy- 
namics is given by the Hamilton equations 



dxo 

It 



Pa ■ 



dpo 
dT 



-w'^xq . 



The general solution can be written as 



xo{T) 
Po{T) 



1 



2uj 



where a and its complex conjugate are fixed by the 
initial conditions. 

Consider next a phase space {x^p) with Hamiltonian 
h — F{hQ), i.e. a function of the standard Hamilto- 
nian for the harmonic oscillator. For instance, the case 
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F{ho) = /iq s-rises in the context of quantum optics, in 
relation with the propagation of hght in non-Unear Kerr 
media . The equations of motion now read 



dx 
dp 



= {x,F{h„)}^F'{ho)p, 

= {p,F{ho)} = -LO^F'{ho)x, 



where F' denotes the derivative of F with respect to its 
argument. These new equations can be easily solved by 
means of a change of time. Specifically, making use of the 
time independence of the Hamiltonian ho (ho — wa^ a on 
solutions to the equations of motion), we can introduce 
the new time parameter T{t) = F'{ho)t. Then, the func- 
tions x{t) = xo[T{t)], p{t) — pi^\T{t)] and the new time 
T serve us to transform the Hamilton equations into the 
standard ones corresponding to the harmonic oscillator. 
So, the classical solutions are 



x{t)^xo[T{t)]-- 
p{t)=po[T{t)]: 



1 



2uj 
iuj 



— iujF' {uja^a)t A iujF' {uja^a)t 

1X0 yJi O 



What we find is an energy dependent redefinition of time 
that induces a different time change for each solution to 
the equations of motion. 

The situation in the quantum theory is quite different. 
The reason lies on the fact than in quantum mechanics 
a physical state does not need to have a definite energy. 
Then, a "change of time" of the form T = F'{ho)t has 
non-trivial consequences for the dynamics. 

The usual quantum theory for the harmonic oscillator 
can be described by introducing a Fock space with cre- 
ation and annihilation operators a^, a. Every initial state 
can be expressed as \(l){0)) — J2'^=q ^nln) , where c„ are 
Fourier coefficients (with the convenient normalization) 
and \n) are energy eigenvectors (that is, ho\n) — nw\n)). 
In the Schrodinger picture, evolving with ho, we find 



\MT)) = e 



-ihnT 



\m) 



However, if the same state \<p{0)) evolves in time accord- 
ing to the evolution generated hy h — F{ho), we get^ 



■iF(nw)t I 



ji=0 



Hence we do not recover an analogous situation to that 
found in the classical system by replacing (formally) the 



^ Notice that the operators ho and h = F{ho) act on the same 
Hilbert space. Moreover, if we demand the function F to have a 
unique absolute minimum at 0, the two Hamiltonians have also 
the same vacuum. 



time T in |0o(T)) with T = F'{ho)t, because |(/'f(0) ^ 
\(j)o{T{t))) except for finear homogeneous functions F. 
Moreover, the properties of the states \(j)o{T)) and \(j)F{t)) 
are quite different. For example, if we consider the 
bounded Hamiltonian h — F{ho) = 1 — e~^° it is ob- 
vious that the high energy contributions to |0f(^)) 
essentially frozen in time with respect to the evolution of 
the low energy ones, in the sense that the phase of the 
former type of contributions remains practically coherent 
in time. 



B. Functionally Related Hamiltonians and 
Cylindrical Waves 

Let us turn now to the discussion of our gravitational 
system. We will deal with the Einstein-Rosen family of 
linearly polarized gravitational waves [17j. These waves 
display the so-called whole cylindrical symmetry [T8| . 
namely, they correspond to topologically trivial space- 
times which possess two linearly independent, commut- 
ing, spacelike, and hypersurface orthogonal Killing vector 
fields. It is well known |^ ^| that these spacetimes ad- 
mit coordinates (T, R, 9, Z) such that the metric is given 
by 



ds' 



{-dT^ + dR^ 



R^de^ 



where ^ and 7 are functions of only R and T. When the 
Einstein equations are imposed, tp encodes the physical 
degrees of freedom and satisfies the usual wave equation 
for an axially symmetric massless scalar field in three- 
dimensions: 

it 

The metric function 7 can be expressed in terms of this 



field on the classical solutions 



expre; 
111 



One gets 



7(i?) 



7o 







1 
2 
1 

2./0 



dRR [[dT^f + {dn^jf] , 
dRR [{drijf + [dR^f] . 



(9) 



Note that 7(i?) and 700 are the energy of the scalar 
field in a ball of radius R and in the whole of the two- 
dimensional flat space, respectively. Furthermore, 700 
coincides with the Hamiltonian Hq given in Eq. {Tj) 2]. 

Nevertheless, to reach a unit asymptotic timelike 
Killing vector field in the actual four-dimensional space- 
time, with respect to which one can truly introduce 
a physical notion of energy (per unit length) 0, [T5| . 
one must make use of a different system of coordinates, 
namely {t,R,9,Z) where T = e~'''°°/^i. In these new 
coordinates the metric has the form 



ds' 



Assuming as a boundary condition that the metric func- 
tion ij] falls off sufficiently fast as i? — > 00, the above 
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metric describes asymptotically flat spacetimes with a, 
generally non-zero, deficit angle. In this asymptotic re- 
gion dt is a unit timelike vector. The Einstein field equa- 
tions can be obtained from a Hamiltonian action princi- 
ple .l^ .15i] where the on-shcU Hamiltonian is given in 
terms of that for the free scalar field by 



H = E{Ho) = 2{l 



.-Ha/2 



(10) 



Owing to these reasons, in the following we will refer to 
t as the physical time and to H as the physical Hamilto- 
nian. 

When we use the T-time and impose regularity at the 
origin i? = , the classical solutions for the field ip can 
be expanded in the form 



m,T)= ^MRk)[A{k)e 



-ikT 



A''{k)( 



ikTl 



A{k) and A^{k) are fixed by the initial conditions and are 
complex conjugate to each other, because ip and Jq (the 
zeroth-order Bessel function of the first kind) are real. 
From Eq. ©, we then obtain 

7oo=Ho= / dkkA'<{k)A{k). 
Jo 

Using this formula, we can express the field in the t-frame 
as 



'dk 



V2 



Jo{Rk) A{k)e 



-ikte 



+A\k} 



ikte-'"=°/^ 



Notice that V'(^,0) = ■(/'ij(^,0). 

In principle, the quantization of the field ^ can be 
carried out in a standard way. We can introduce a 
Fock space in which ip(R,0), the quantum counterpart 
of ^{R,0), is an operator-valued distribution Its 
action is determined by those of A{k) and A'''{k), the 
usual annihilation and creation operators, whose only 
non-vanishing commutators are 



i(fcl),it(fc2) =J(fci,fc2) 



Explicitly, 



°dk 



^P{R,0)^4'EiR,0)^ -^MRk) A{k) + AHk) 
Jo v2 L 



(11) 



(12) 



In the Schrodinger picture, operators do not evolve, and 
the problem of the time evolution is transferred to the 
physical states. We will return later to this issue. In the 
Heisenberg picture, on the contrary, operators change in 
time and states remain fixed. In this case, the value of 
the quantum field at any time can be obtained from 
its value at T = t = 0, evolving in one of the two times 
that we have at hand. One is the physical time t, with 
associated Hamiltonian H given in Eq. H1U|) . The other is 



the time T of the auxiliary three-dimensional Minkowski 
spacetime, its Hamiltonian being that of a massless scalar 
field, Hq. From our discussion in Sec. II, this time can 
be identified with the perturbative time (denoted also by 
T) that arises in the study of the Einstein-Rosen waves 
in linearized gravity. 

In the perturbative time T, the evolution is provided 
by the unitary operator Uq{T) — exp(— iTifo) where 



= / dkkA\k)A{k) 

'0 



(13) 



is the quantum Hamiltonian of an axially symmetric 
scalar field in three dimensions. Then 



V'(i?,r)= c/j(r)«'(i?,o)c/o(r) 

'dk 



V2' 



MRk) A{k)e-'''^+A^k)( 



ikT 



(14) 



The quantization procedure in this case is very simple: 
we have substituted the initial conditions A{k) and A''{k) 
in the classical solution by the corresponding quantum 
operators. 

The situation changes when we choose the physical 
time t as the time parameter. The quantum Hamilto- 
nian can be defined as H = E{Hq) = 2(1 - e^"°^^). 
We can then reach a unitary evolution by means of 
U{t) = exp{—itH). This leads to the following time 
evolved operators 



AE{k,t) := U^{t)A{k)U{t) = cxp\~itE{k)e'""^'^ 



iyfc,t) = A\k) exp\itE{k)e-"°^^ 



A{k), 
(15) 



where E{k) — 2(1 — e^'^^^). It is important to realize 
that the quantum evolution in the physical i-frame is 
not obtained by changing A{k), A''{k), and 700 by their 
direct quantum counterparts. In fact, by restoring the 
value of the dimensionful constants h and Ga, we can 



write E{k) = (1 



-4Gfe 



)/(4G3) with G = hG3, so that 



-E(k) = tk + o{h), 

h 



and we can expand ^^(fc, t) and A*^{k, t) in powers of h, 
AE{k,t) = exp(-itke-'^'^"A A{k) +o{h), 



Al{k,t) ^ A\k) cxp Utke-'^^^") + o{n) 



(Here, we have expressed Hq with dimensions of an in- 
verse length) 

Setting again 8G — 1, the expansion above clearly 
shows that the quantum evolution of the creation and 
annihilation variables in the physical time diff'ers from 
the "classical evolution" in higher-order quantum correc- 
tions. 

This unusual behavior can be partially corrected in the 
sense that one can actually find an operator for H such 
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that the quantum evolution in the t-time is similar to the 
classical one. In fact, if one considers the normal ordered 
Hamiltonian Hnor —■ H : and its associated unitary evo- 
lution operator Unor — exp(— iii7„or), it is easy to prove 
that 

Anor{k,t) ^ exp (^-itk : e^"°^^ A{k), 
Aiorik,t) - it (A:) exp (^itk : e^"/^ . 

Therefore, the quantum evolution in t given by Hnor par- 
allels that encountered in the classical picture. Unfor- 
tunately, there is a severe problem with this quantum 
dynamics that renders it physically unacceptable: the 
Hamiltonian Hnor is unbounded both from above and 
below 2]. 

Obviously, the different possibilities for the unitary 
time evolution of the field ip considered here are closely 
related. Specifically, there exists a unitary mapping be- 
tween both types of evolution. This is a consequence of 
the fact that they are both unitary and coincide with 
the identity in the same Hilbert space at T = i = 0. 
So, if we denote the operators evolved from X(0) in 
the times T and t by X{T) = ul{T)X{0)UoiT) and 
Xsit) = UHt)X{0)U{t), respectively, we obtain 

XEit) = U\0o{T)X{T)Ul{T)U{t) . 

Then, we can go from one evolution to the other by means 
of the unitary operator Ul{T)U{t). 

Let us close this section with a few comments about 
the Fock space on which the field i> acts. In the Heisen- 
berg picture the states do not depend on time. They 
are constructed by successive actions of the time inde- 
pendent creation operators on the vacuum of the theory. 
It is important to point out that both Ho and H act on 
the same Hilbert space and have the same vacuum, that 
will be referred to as |0). Explicitly, given any square in- 
tegrable complex function 0„(fci, . . . , fc„) (with the con- 
venient normalization) we can write a n-particle state in 
the form 

\(pn)^ dkr ■ ■ / (fci, . . . , kn)A^{ki)- ■ ■ it(fc„)|0). 

Jo Jo 

According to the usual interpretation of quantum me- 
chanics, the measurable physical quantities correspond 
to expectation values of observables. We can go over to 
the Schrodinger picture by assigning the time evolution 
to the quantum states 

x{t; 0) = {m\u\t)x{o)u{tMo)) Mmixmm- 

Defining \Mt)) - ^>WI0(O)), |<^(T)> - Uo{T)\m), 
and noticing that the U -operators satisfy 

idtU{t) = HU{t), idrUoiT) = HoUo{T) , 

it is straightforward to see that the evolved states are 
solutions to the Schrodinger equations 



As in the Heisenberg picture, the unitary operator 
lJ{t)ijl{T) provides the bridge between the two kinds 
of quantum evolution, 

\cpEit))=U{t)ul{TmT))- 

Particularizing the above results to the case of n-particle 
states we readily get that, for the T-time, 

|0„(T)> -J7o(r)|0„) = 

/ dkr ■ ■ dfc„e-^'=""^0„(fci, kn)A\ki)- ■ ■ A\kn)\0) 
Jo Jo 

where ktot = X]J=i ^j- Notice that |0„(r)) is a superpo- 
sition of eigenvectors of Hq with eigenvalues equal to ktot- 
On the other hand, if we evolve the states with U{t), it 
is not difficult to check that 

/ dkr ■ ■ dfc„e-^^('='-)*0„(fci, kn)A\kiy ■ A\knm. 
Jo Jo 

In other words, 10b. n(i)) is a superposition of eigenvec- 
tors of H, each of them with energy E{ktot)- Finally, it 
is worth pointing out that the iJ-energy is not additive: 

n n 

E{ktot) = 2- 2 e-(^? VE^(^*) = 2^ - 2^e-^'/2 

i=l i=l 

This property is directly related to the existence of an 
upper bound for the physical Hamiltonian. 

IV. MICROCAUSALITY 
A. Free Hamiltonian 

Microcausality plays a crucial role in perturbative 
quantum field theory; in fact it is a crucial ingredient in 
such important issues as the spin-statistics theorem jl3j]. 
The point of view that we will develop in this section is 
the idea that field commutators of the scalar field that en- 
codes the physical information in linearly polarized cylin- 
drical waves can be used as an alternative to the metric 
operator to extract physical information about quantum 
spacetime. Some relevant information may be lost, but 
the availability of explicit exact expressions for commuta- 
tors [even under the evolution given by the Hamiltonian 
(|1U|| ] opens up the possibility of getting precise informa- 
tion about the quantum causal structure of spacetime. In 
particular we will see how the light cones get smeared by 
quantum corrections in a precise and quantitative way. 

As is well known, the study of causality in perturbative 
quantum field theory requires the consideration of mea- 
surements of observables at different spacetime points 
and their mutual influence. The key question is whether 
measurements taken at spatial separations commute or 
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not. The relevant commutators can be seen to be pro- 
portional to those of the quantum fields at two different 
spacetime points x and y. For a free massless scalar field 
described by the standard Lagrangian, this commutator 
is a (c-number) function of x and y that is exactly zero 
when x — yis space-like [l^. This means that observables 
at points separated by spatial intervals commute. Some- 
thing similar happens for fermions described by the usual 
Lagrangians if, instead of commutators, one takes anti- 
commutators of the fields (observables for fermion fields 
can be written as even powers of them and the relevant 
commutators can be written in terms of anticommutators 
of the basic fields) |12f |. 

In the case that we are considering in this work the 
only physical local degree of freedom that we have is given 
by the scalar field V'- What we will do in the following 
is discussing microcausality by considering the different 
time evolutions introduced in the first part of the paper. 

We start by computing [■(/j(_Ri, Ti), '0(i?2, TIj)] for the 
field operators 1)14(1 obtained in the Heisenberg picture 
by evolving the field at time T — Q with the Hamiltonian 
Hg given in Eq. l(T^ . It is straightforward to get (see 
[20| for a somewhat related computation) 

[i^(i?i,ri),v>(i?2,r2)] = 

i / dkJoiRik)Jo{R2k)sm[{T2-Ti)k]. (16) 
Jo 

Let us discuss the main features of this commutation 
function, that we will refer to in the following as the 
iJo-commutator. To begin with we can easily see that, 
as it happens for the familiar free field theories, the com- 
mutator H16|) is a c-number, i.e. it is proportional to the 
identity in the Fock space. In addition, for Ri fixed, let 
us call regions I, II, and III the regions of the {R2, T2—T1) 
plane defined, respectively, by < IT2 — Ti| < |i?2 ~ -RiL 
\R2-Ri\ < \T2-Ti\ < R2 + R1, and R1+R2 < \T2-Ti\. 
Then, it can be shown that the i/o-commutator vanishes 
in region I, whereas in region II it can be written as pH 

[^(i?i,ri),v-(i?2,r2)] = 

i 1 



T2-T, 



7tVWR2 



K 



{T2-t^Y-{R2-RiY 
ARIR2 



(17) 



Finally, its value in region III is |21 



[V;(i?i,Ti),^(i?2,T2)] 
2i 



K 



4i?ii?2 



V'7r2[(T2-Ti)2-(i?2-i?i)2] \^{T2-T,f-{R2-Ri? )' 

Here, K denotes the complete elliptic integral of the first 
kind, K{k) := /J'^^dg/Vl - P gin^ 6* (alternatively, the 
integral in Eq. H16|l can be written in terms of the asso- 
ciated Legendre functions P_i/2 and Q-1/2 HJl)- 

Several plots of the absolute value of this function (in 
fact, of its imaginary part) for fixed values of Ri and 




R2 



FIG. 1: Regions in the (-R2, 72 — Ti) plane for the ffo-commu- 
tator. Region I corresponds to < jr2— Ti] < \R2 — Ri\, 
region II to |i?2 - < lr2 - Ti] < J?2 -I- -Ri, and region III 
to R2 + R\ < |r2 — The singularity of the Ho-commu- 
tation function lies in the boundary between regions II and III, 
whereas the singularity for the _E(ffo)-commutator appears 
for Ri= R2. 





- T2-T, 




T2-T1 



T2-T1 



FIG. 2: Absolute value A of the commutator of the scalar 
field [V)(-Ri, Ti), i/;(i?2, Tb)] for R\ = 1 and several values of 
i?2 as a function of T2 — Ti . 



i?2 are shown in Fig. |21 and a three-dimensional plot 
for a fixed value of i?i can be found in Fig. O Some 
interesting properties can be read off from the previous 
expressions. First of all we see that the commutator is 
in fact identically zero in the region labeled I; outside it 
differs from zero. Second we see that it is singular in 
the line T2 — Ti = R2 + Ri- It is easy to check that the 
singularity is logarithmic by using the explicit form of the 
propagator (|17|) and expanding around T2—T1 = i?2-l--Ri. 
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FIG. 3: Absolute value A of the commutator of the scalar field 
[i){Ri , Ti), V;(7?2, 12)] for i?i = 1 as a function of T2 - Ti and 
R2 ■ The light cone structure and the singularity at T2 — T1 = 
R2 + Ri can be readily seen. The plot for negative values of 
T2 — Ti is obtained by reflecting with respect to the R2 axis. 



This singularity can be understood as a consequence of 
the presence of the symmetry axis; this is supported by 
the fact that the _ffo-commutation function satisfies 



d-^-^dR{RdR)]-mR,T),^{R2,T2)]=0, 



R 



dniRdR)]- e{T -T2)[iP{R,T),i,{R2,T2 
-^S^^\R- R2,T~T2), 



which shows that the full commutator is an axially sym- 
metric solution to the 2+1 dimensional wave equation 
and the future part of it (obtained by multiplying by a 
step function) is an axially symmetric Green function for 
the same equation. It is important to point out that 
even though there is a Minkowskian background metric, 
the presence of a center of symmetry breaks Lorentz in- 
variance; in fact, the only symmetries of the 2+1 axially 
symmetric wave equation 



^dR{RdR)j^T,R) = 



that transform solutions <f>(r, R) into solutions'^ are 
i) Translations 

$(T, i?) $(T - e, i?) (eeK), 



ii) Dilatations 

<i>{T,R)^ <^>{e'T,e'R) (e e M), 

iii) Inversions 
<I>(T, R) ^ '^[r.(T,R),p.iT,R)] ^^^^^ ^^^^^^^ 



^i-2re+(r2-i?2)£2 



with 



TeiT,R) = 
Pe{T,R) = 



T~e{T^-R^) 
l-2Te+ (r2 - i?2)e2 
R 

l~2Te+ (r2 - i?2)£2 



These can be obtained systematically by using the gen- 
eral theory of symmetry groups for partial differential 
equations (see p2l|V 



B. E{Ho) Hamiltonian 

Let us consider now the commutator of the field oper- 
ators obtained in the Heisenberg picture with the quan- 
tum Hamiltonian H E{Hq) = (1 - e-'^^^o)/{4:G3). 
Here G = G^H, and we have restored the values of the 
Planck constant and the three-dimensional gravitational 
constant (but kept c — 1) to have the possibility of dis- 
cussing the semiclassical limit. Recall that the Hamilto- 
nian Hq is given by Eq. (|13|l . where we choose the nor- 
malization of the creation and annihilation operators so 
that the commutation relation (|ll|l is satisfied (with no 
h in the r.h.s.). Notice that, with this convention, A{k), 
Hq, and H have formally the same dimensions as (5^/2 (or 
A:~^/2)^ (7-1 (or hk), respectively. To distinguish 

the field operators evolved with the new Hamiltonian H, 
we will denote them by ipEiRTt). In order to compute 
the commutator [■0_e(^1j ^£;(^2j i2)]j we make use of 
the expressions (|15f) for the creation and annihilation op- 
erators evolved in t, which can also be written as 



Asik^t) = A{k)exp \it l^E{Ho - k) - E{Ho) 
iyfc,t) = it(A:)exp it{E{Ha + k)- E{Ho) 
Substituting then the relation 

/"OC 

i^E{R,t) ^ Vid dkJo{Rk) AE{k,t) + A^j^{k,t) 
Jq l 

we get the relevant field commutator 

/>oo />oo 

[i^E (Rl,tl)jE{R2,t2)]= dkJ dk2 MRl kl)Jo {R2 k2 ) 

Jo Jo 

x4G\[AE{ki,h),AEik2,t2)] + [A^Eiki,ti),A^Eik2,t2)] 



(18) 



^ In addition to the ones coming from its linear and homogeneous 
character. 



-[AE{ki,h),Al{k2,t2)] + [Al{ki,ti),AE{k2,t2)] 
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G = 0.02, t2-ti = 0.6 



G = 0.011, t2-ti = 0.6 




R2 




R2 



G = 0.02, t2-ti = l 



G = 0.011, t2-ti = l 




FIG. 4; Comparison between the absolute values (over 8G) 
of the i/o-commutator and the vacuum expectation value of 
the i?(iifo)-commutator for two different values of G, plotted 
for Ri = l. 



The commutators involving AE{k,t) and can 
be found in Appendix ^ Note that, in contrast with 
the evolution given by Hq, for which the commutation 
function is a c-number, the situation now is more com- 
plicated because [ilJE{Ri,ti),ipE{R2,t2)] is an operator, 
as it happens in interacting theories. We will analyze two 
types of matrix elements for it, namely, the expectation 
value on the vacuum and on one-particle states. In ad- 
dition, we will briefly comment on the expectation value 
on coherent states. 



Vacuum expectation value. 



[0\[^E{Rl,h),ijE{R2,t2)]\0) 

iSG dkJo{Rik)Jo{R2k)sm 
Jo 



t2-h 

4G 



(l-e-4Gfc) 



(19) 



In the following we will discuss the differences and simi- 
larities of this matrix clement and the iJo-commutation 
function. 

We point out that the factor sin[(r2 — Ti)fc] that ap- 
pears in the integrand of the iJo-commutator is substi- 
tuted now by sin[(t2 - ti)(l - e-^^^) / {AG)]. They coin- 
cide for fc — > 0, but the former of these functions oscillates 
for all values of k whereas the latter approaches a con- 
stant value when fc — > oo. This changes the convergence 
properties of the integral. In particular it is straight- 
forward to see that the integral (|19|l diverges whenever 
Ri = i?2 (except if sin [(^2 - ti)/{4:G)] = 0) but con- 
verges otherwise. Therefore, the vacuum expectation 
value has a singularity structure that differs from the one 



given by Eq. H16(l . This has some interesting physical 
consequences. First we see that the singularity that orig- 
inates in the axis in the Hq case is not present when the 
evolution is dictated by E{Ho); this can be interpreted as 
a blurring of the axis due to quantum corrections. Sec- 
ond we see that a completely different kind of singularity 
pops up when the evolution is generated by E{Hq). From 
a mathematical point of view its origin is clearly related 
to the fact that the energy is bounded from above and, 
hence, for large values of k the integrand is just a product 
of two Jo Bessel functions (which give a divergent integral 
if their arguments coincide). Physically, the emergence 
of the singularity can be understood in an intuitive man- 
ner by writing a state as a superposition of vectors of the 
form i^(fc)|0), 



i^E{RiMM^- 



; / dfc Jo(i?ifc)e^(i-"""''Ut(fc)|0), 
^0 



and projecting onto V;i=;(i?2,t2)|0). The e'*''-'^-" "'')/(4G) 
factor (for j = 1 or 2) goes to a phase that depends only 
on tj as fc — > oo. So, if Ri = R2, the coefficients of the 
linear superposition defining, respectively, 'ipE{R,ti)\0) 
and ^I^e{R, ^2)10) differ only by a constant phase for large 
values of fc. This means that, in the sector of large fc, 
these two states have coherent phases and therefore a 
constructive interference. Since each of them has an in- 
finite norm, their scalar product diverges. A similar ef- 
fect can be found in the quantum dynamics of a free 
particle with an energy given by a function E{k) with 
dE{k)/dk — > as fc cx). If one builds a wave packet as 



dkC{k)e 



i[tE{k)+kx 



with C(fc) peaked around a large value of fc, fco 3> 1, 
the group velocity becomes almost zero and $(<:, x) stays 
essentially the same at every x for long periods of time. 

The integral in Eq. H19|) can be written explicitly as a 
convergent series (when Ri ^ R2) by expanding the sine 
function as a power series of e~'^^^ and computing the 
resulting integrals involving two Bessel functions and an 
exponential 21] 



['4^E{RlM),i^E{R2,t2W) = 

(-l)"(At)2 



i8G 



;sm 



(A^)E 



n=0 



(2n)! 



■Q-1/2 ['^2n{Rl, R2)] 



n/ A -i\2n+l 



- cos(At)E ^ ]2n+%. i?2)] !> . (20) 



n=0 



{2n + 1) 



Here, we have introduced the notation 
A. = ^, .„(i?„i?2) = i'^'^' + ^? + ^^' 



4G 



2R1R2 



Besides, Q_i/2(a;) = ttF (|, i; 1; ^) /v^ [with x > 1] 
is the associated Legendre function of the second kind 
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[2l|. We recall that the function Q-i/2{x) grows without 
bound as the argument approaches x = 1 and falls off 
to zero as 7r/V2x when x ^ oo. When Ri = R2 the 
singularity in Eq. (|20|l comes just from the n = term 
in the first series of the expansion, that is given by 



Ri = i?2 singularity. This can be seen by considering the 
last term in Eq. H21I) : the integral in ^2 is 



iSG 



n,/RjT.. 



Q- 



1/2 



Ri + R2 
2R1R2 



t2-h 

4G 



A series of plots of [Vi(i?i,Ti = ti), ■0(i?2, 12 = ^2)] and 
(0| [V'£;(i?i,ti),V'ij(i?2,t2)] |0) (both over 8G) is shown in 
Fig. ^ for fixed values of Ri and t2 — ti as a function 
of i?2, with several choices for G. We choose t2 — ti 
small enough to guarantee the rapid convergence of the 
series in Eq. (|2U|) and leave a discussion of the behav- 
ior of the integral H19I) when G ^ for future work. 
As we can see (0| [?/;_e(-Ri, ii), '0£;(^2, ^2)] |0) seems to ap- 
proach [ip{Ri,Ti — ti),ip{R2,T2 — t2)] at least in a cer- 
tain average sense when G is sufficiently small (though 
not vanishing). It falls off to zero quite quickly outside 
the light cone defined by the free commutator and the 
auxiliary Minkowski metric, and displays an oscillatory 
behavior within this light cone. The characteristic length 
of this oscillation decreases with G, as well as close to 
the i?i — i?2 singularity. The approximation obtained 
by truncating the series expansion H20|l . keeping a suf- 
ficiently large number of terms, compares well with the 
results of numerically computing expression H19() , at least 
for low enough values of t2 — ti. 

• Expectation value on one-particle states. 

We consider now states of the form 

\p) = / dkf{k)A\k)\0) 
Jo 

where the function / = \f\e^'^f satisfies J^dk\f{k)\'^ — 1. 
We then have 

{p\\i,E{RiM),i>E{R2M)\\p) - 

~i8G Jdki Jdk2MRiki)hjo{R2k2)\.f(k,)f{k2)\ 

X sin [2G{ti - t2)E{ki)E{k2)] cos [17/(fci, fca)] 

- Jo(i?2fci)|/(fc2)|' sin [(^2 - ii)£;(fci)e-4'5^^1 1, (21) 



where 

Q.f{ki,k2) 



2G{t2 - h)E{ki)E{k2) + tiE{ki) 

t2E{k2)~(bf{ki)+(f>f{k2). 



A complete discussion of the meaning of the previous 
expression is beyond the scope of this paper. Neverthe- 
less, some features already present in the vacuum ex- 
pectation value are also present here; in particular the 



dfc2|/(fc2)|'sin {t2-h)E{ki)e 



-4Gfc2 



which takes in general a non-vanishing constant value 
(depending on t2 — ti and G) as fci ^00, thus rendering 
the remaining integral in ki divergent. As the first term 
in Eq. (|21() leads to a convergent integral, we conclude 
that the expectation value is singular when Ri = i?2- 

It is not difficult to obtain as well an explicit expres- 
sion for the expectation value of the ii^(7?o)-commutator 
on the coherent states of the field ip. These diagonal 
matrix elements are calculated in Appendix B. For our 
discussion in this work, let us only comment that the re- 
sult is actually divergent when Ri = i?2- This supplies 
further support to the claim that the considered singu- 
larity is indeed a generic feature of the system. 



V. CONCLUSIONS AND PERSPECTIVES 

Linearly polarized cylindrical waves can be studied in 
great detail both from the classical and quantum points 
of view. As we have seen, there are two relevant Hamil- 
tonians for the study of the system. We have shown that 
the action and the metric of the gauge-fixed model in 
linearized gravity reproduce the results obtained by con- 
sidering full cylindrical gravity and working to the first 
perturbative order. We get in this way a free Hamilto- 
nian. The Hamiltonian governing the dynamics of the 
full system, on the other hand, is different from the free 
one, but turns out to be a function of it and presents cer- 
tain features with deep physical consequences, like e.g. 
the existence of an upper bound. 

We have studied the similarities and differences of 
these two admissible kinds of evolution; in particular, 
we have discussed how the emergence of an upper bound 
for the energy affects the causal structure of the model 
and the spreading of the light cones. The field commu- 
tator for the free Hamiltonian is a c-number and shows 
the typical light cone structure found in standard per- 
turbative quantum field theories. The commutator for 
the physical Hamiltonian, as it usually happens for in- 
teracting theories, is no longer a c-number, so one has 
to consider its matrix elements. By concentrating on the 
vacuum expectation value we have been able to see sev- 
eral interesting phenomena: a spreading of the light cone 
as a function of the gravitational constant, the disappear- 
ance of the singularity present in the free case due to the 
smearing of the symmetry axis and the appearance of a 
new type of singularity associated with the fact that the 
energy is bounded from above. This new singularity is 
also present in the other expectation values discussed in 
the paper, namely for one-particle states and coherent 
states, and appears to be a generic feature of the model. 

There are several open questions that we plan to ad- 
dress in future work. In particular, it would be desirable 
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to reach a better understanding of the behavior of the 
field commutator in the hmit in which the length scale 
provided by G goes to zero. The expectation values of 
the i?(_ffo)-commutator discussed here resemble those de- 
rived from the free Hamiltonian Hq at least in a certain 
average sense. However it is not obvious how precisely 
and up to what extent they actually relate to each other. 
This is partly so because of the different singularity struc- 
ture found in both cases. Further research on this sub- 
ject will concentrate on the properties of the model in the 
semiclassical limit G ^ 0. We will also pay detailed at- 
tention to matrix elements of the field commutator other 
than the vacuum expectation value, with the aim at dis- 
cussing how the smearing of the light cones depends on 



the energy. 



APPENDIX A: USEFUL COMMUTATORS 

In this appendix, we compute the commutators of 
the creation and annihilation operators AE{k^t) and 
obtained from the corresponding operators 
A{k) and {k) via the unitary evolution generated by 
E{Ha), where Hq = dkkA^k)A{k). Employing rela- 
tions H18|l and the basic commutators (|11() . it is possible 
to show that 



[AEikuti),AE{k2,t2)] 
[iyfcl,tl),iyfc2,t2)] 
[AE{ki,h),AUk2,t2)] 

[Aliki,ti),AEik2,t2)] 



A{ki)A{k2) exp itiE{Ho - fci - fca) + i(t2 - h)E{HQ - /ca) - it2E{Ho 
A{ki)A{k2) exp it2E{Ho - fci - ^2) + i{ti - t2)E{HQ - ki) - itiE{HQ 
it(fci)it(fc2) exp \ihE{Ho + ki + /ca) + i{t2 - ti)E{Ho + fca) - it2E{Ho) 
A^k,)A\k2) exp [it2E{Hn + k,+ k2) + i{h - t2)E{Ho + fci) - itiE{Ho) 
itiE{Ho + k2- fci) + iit2 - h)E{Ho + fc2) - it2E{Ha) 
it2E{Ho + k2- fci) + i{h - t2)E{Ho - fci) - thE{Ho) 
itiE{Ho + fci - fca) + i{t2 - ti)E{Ho - fc2) - i<2£^(i?o) 
it2E{Ho 4- fci - fc2) + iih - t2)E{Ho + ki) - itiE{Ha) 
I 



i(fci)it(fc2)exp 
it(A:2)i(fci)exp 
it(fci)i(fc2)exp 
i(fc2)i^(fci)exp 



APPENDIX B: EXPECTATION VALUES ON 
COHERENT STATES 

We consider coherent states of the field ip, given by 

dk 



\^c) = if c exp 



/8G 



C(fc)At(fc) |0) 



= -c|l(f^C(.)..,.,)-|„,, 

where C(fc) is a square integrable function and Kc is a 
normalization constant satisfying 



where 

/.(G)=exp|^°'g|C(fc)p(e-4S^'=-l) 

G,(fci, fc2) = 8G^(fci, fc2)|[6(fcl, -fc2)] -[6(fc2, -fcl)]^ 

+C(fc2)G(fci)| [5(fci, -fc2)]^- [6(-fc2, fci)]^ |e-4G^'=i 
+G(fci)G(fc2)|[fc(-fci,fc2)]^-[6(fc2,-fci)]^je-^^^'=^ 
+G(fci)G(fc2)l [Hfci, ^2)]^ - [fo(fc2, fci)]4e"''^'^''^''^ 



|ifcP = exp(^-^°°g|G(fc)|^ 
The expectation value of the i?(_ffo)-commutator is 

{^c\[4'EiRl,h),4>EiR2,t2)]\'^c) = 

/ ^ / dfc2 Jo fci ) Jo (i?2 fc2 ) V ^ G, (fci , fc2 ) , 
^0 ^ ^0 s- 



-fG(fci)G(fc2) [fo(-fcl,-fc2)]-^-[fo(-fc2,-fcl)]' 



and we have employed the notation 



b{kn,km) — 



i„e4«'=-(e4«'="~l) + i„(e^^'^--l) 



4Gfc„ 
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Note that, when i?i = i?2, the dcha in the expression of 
Gs{ki,k2) leads to the divergent integral 



/' 

Jo 



dkJl{Rik) 



8GL 



2s+l 



^(2. + l)! 



^(^2-^l)(l-e-"^'') 
4G 



-i2s+l 
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